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Abstract 

^y. . We introduce discrete time Markov chains that preserve uniform mea- 

sures on boxed plane partitions. Elementary Markov steps change the size 
of the box from a x b x c to (a — 1) x (b + 1) x c or (a + 1) x (b — 1) x c. 
Algorithmic realization of each step involves 0((a + b)c) operations. One 

\^ ' application is an efficient perfect random sampling algorithm for uniformly 

f) , distributed boxed plane partitions. 

Trajectories of our Markov chains can be viewed as random point 
configurations in the three-dimensional lattice. We compute the bulk 
C^ , limits of the correlation functions of the resulting random point process 

on suitable two-dimensional sections. The limiting correlation functions 
define a two-dimensional determinantal point processes with certain Gibbs 
properties. 
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Intoduction 



For any integers a,b,c> 1 consider a hexagon with sides a, b, c, a, b, c drawn on a 
regular triangular lattice. Denote by fl a xbxc the set of all tilings of this hexagon 
by rhombi obtained by gluing two of the neighboring elementary triangles to- 
gether (such rhombi are called lozenges). Equivalently, fl a xbxc is the set of all 
dimers on the part of the dual hexagonal lattice cut out by our (a, b, c)-hexagon. 
An element of H4 x 5 x 5 is shown on Figure 1 . 

Elements of fl a xbxc have a number of different interpretations, see e.g. Sec- 
^^ ' tion 1 below. In particular, they can be viewed as plane partitions or stepped 

surfaces inside a three-dimensional box of size axbxc. A uniformly distributed 
element of Qaxbxc then provides a basic model of a random surface. This 
model and its generalizations have been thoroughly studied, see e.g. [CLP , 
[CKP] . [DMB] . [Ei], fjl] [J2], [JN], Hil], [KO], [EE], [EL], [P2], [Wl], [W2]. 
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Figure 1. Tiling, stepped surface or boxed plane partition. Separating lines 
between horizontal lozenges are removed. 

The main goal of this paper is to introduce and study certain discrete time 
Markov chains on boxed plane partitions that preserve the uniform measures. 

Denote by fi a xbxc the uniform probability measure on £l a xbxc- We construct 
two families of stochastic matrices 



P 



axbxc 



a 



xfcx 



c X fya-1 



(a-l)x(6+l)xc 



[0,1] 



such that 



Paxbxc ' ^axbxc X f^(a+l)x(6-l)xc ~~ * [0, 1] 



22 Max6xc(w)P ax6xc (cj,w') =M(oTl)x(6±l)Xc(w') 



wef2„ X i, Xc 

for all a, 6, c such that the participating sets fi are nonempty. 

Although it is a little awkward to write matrices P in one formula, ap- 
plication of Markov operators corresponding to these matrices is fairly easy 
to describe algorithmically. The exact algorithm can be found in Section SJ 
Roughly speaking, the algorithm does the following: Given uj G Q a xbxc, in or- 
der to construct u/ distributed according to -P OX (, Xc (<*>,u/) it needs to consider 
all horizontal lozenges of u> sequentially from left to right. For each such lozenge 
the algorithm decides on its new position using a simple one-dimensional prob- 
ability distribution. (Here we ignored occasional appearance/diappearance of 
horizontal lozenges on top and at the bottom, that sometimes take place.) In a 
way, this means that P^ decompose into products of one-dimensional Markov 
steps. 

The algorithm has some similarity to the shuffling algorithm for domino 
tilings of the Aztec diamonds introduced in EKLPJ. Indeed, that algorithm 
also maps uniform measures to uniform measures (actually, it works for a one- 
dimensional family of measures that includes the uniform one), and it also de- 
composes into simple (Bernoulli) Markov steps. For that reason we call our 
algorithm the shuffling algorithm for boxed plane partitions. 



We further consider Markov chains obtained by successive application of 
arbitrary sequences of matrices P + and P~ . The initial condition and the 
one-time distributions are all uniform measures on the appropriate spaces Q. 
One example is the application of the sequence of b matrices P + to the unique 
probability measure on singleton ^( a +h) x oxc- This gives a perfect sampling 
algorithm for [i a xbxc that takes 0((a + b)bc) one-dimensional steps. When a, 
b, and c are comparable, the algorithm is roughly as efficient as that from |Kr| , 
and it is more efficient than other known algorithms, cf. |BMj , |LRS] , [Plj . [P2j . 
|Wlj . |W2] . Another example is an alternating sequence (P + P~)(P + P~) ■ ■ ■ 
that provides an equilibrium dynamics on £l a xbxc with [i a xbxc as the equilibrium 
measure. 

These Markov chains can be viewed as two-dimensional random 
growth/decay models. One important feature of these models is the fact that 
on certain two-dimensional sections of the three-dimensional space-time, their 
suitably defined correlation functions are computable in a closed dctcrminantal 
form. For the static correlation functions (those of the measures Haxbxc) this 
fact is well known, see [Kaj , |Ke2] , [J2] , [Gor] . 

We then focus on the bulk asymptotics of the correlation functions. Namely, 
when a,b,c — > oo with a : b : c fixed, we look at the three-dimensional lattice 
process near a fixed point of the global limit shape of our random surface. On 
the same two-dimensional sections we compute the limiting correlation functions 
that also have the determinantal form. In fact, for every such section, they 
define a two-dimensional determinantal random point process with a certain 
Gibbs property, see |BSj and references therein. Our results extend similar 
results for the static case obtained in |Gorj . 

The construction of P is an application of the general algebraic formalism 
developed in |BFj . The continuous time Markov chain considered in |BF] in 
detail can be viewed as the degeneration of P ± near a corner of the hexagon 
as a, b, c become large, and either a and b is substantially larger than the other 
two. It is worth noting that the shuffling algorithm for domino tilings of the 
Aztec diamonds also fits into the formalism of }BF] , see Section 2.6 of |BF| and 

The proof of the bulk asymptotics involves spectral decomposition of the 
matrices P in terms of Hahn classical orthogonal polynomials. However, the 
limiting argument does not require the asymptotics of Hahn polynomials them- 
selves — it only requires the much simpler asymptotics of the difference op- 
erators related to Hahn polynomials. This approach to bulk asymptotics of 
determinantal point processes is due to G. Olshanski; it was first used in [BQ2| 
for Charlier and Krawtchouk orthogonal polynomials, and in [Gor] it was further 
developed in the more complex case of Hahn polynomials. 

The paper is organized as follows. In Section 1 we discuss different com- 
binatorial definitions of our model and introduce notations. In Section 2, we 
introduce four stochastic matrices on one-dimensional point configurations and 
prove certain commutativity relations between them. In Section 3, we use these 
four matrices to define P^; the construction is based on an idea from [DFj . 
Section 4 contains the algorithmic description of P and images obtained from 



their computer realizations. In Section 5 we compute the correlation functions 
on suitable two-dimensional sections of the space-time, and in Section 6 we 
obtain the bulk limits. 

The first named author (A. B.) was partially supported by the NSF grant 
DMS-0707163. The second named author (V. G.) was partially supported by 
RFBR grant 07-01-91209, the Moebius Contest Foundation for Young Scientists 
and Leonhard Euler's Fund of Russian Mathematics Support. 

1 Basic model 

The main object of our study has many different combinatorial interpretations. 
In this section we discuss some of them. 

Consider a tiling of an equi-angular hexagon of side lengths a, b, c, a, b, c 
by rhombi with angles 7r/3 and 27r/3 and side lengths 1. Such rhombi are called 
lozenges. 

Lozenge tilings of a hexagon can be identified with 3d Young diagrams 
(equivalently, boxed plane partitions) or with stepped surfaces. The bijection is 
best described pictorially. Examine Figure 1, where a tiling of the (4x5x5) 
hexagon is shown, and view this picture as a 3d shape. 

Given a tiling we construct a family of non-intersecting paths on the surface 
of the corresponding 3d Young diagram as in LRS] , [J2] . Figure 2 provides an 
example. 




Figure 2. Non-intersecting paths on the surface of 3d Young diagram. 

We view this family as a family of paths on the plane. It is convenient for 
us to do one more modification. We replace the downgoing segments of paths 
by horizontal ones and we replace upgoing segments by segments of slope 1. 
Consequently, our family is interpreted as a family of non-intersecting paths on 
Z 2 . Figure 3 shows the family corresponding to the tiling on Figures 1 and 2. 




Figure 3. Non-intersecting paths on Z 2 . 

Below we are going to use the last interpretation. 

Let us introduce some notations. Denote by N the number of paths (in our 
example N = 5). Introduce coordinates on Z 2 so that the first path starts at the 
point (0,0) and ends at the point (T,S). The second path starts at the point 
(0, 1) and ends at the point (T, 5* + 1), and so on. Finally, the Nth path starts 
at the point (0, N — 1) and ends at the point (T, S + N — 1). In our example 
T = 9 and S = 5. 

Denote by Cl(N,T, S) the set of families of N non- intersecting paths made 
of segments of slopes 0,1, starting at (0,0), . . . , (0, TV — 1) and ending at 
{T,S),...,(T,S + N-1). Note that ft (TV, T,S) ^ 0, is equivalent to < S <T. 

Denote by /i(iV, T, S) the uniform measure on il(N, T, S). 

Note that £l(N,T, S) was called Q a xbxc in Introduction, while measure 
H(N, T, S) was called fiaxbxc- Here a = T - S, b = S, c = N. 

Set 



r^ T = {ieZ: max(0, t + S - T) < x < mm(t + N - 1, S + N - 1)} 



and 



~S,t 
jS.t 



X 



S,t 
N,T 



{{X1,X1, 



, X N ) e (£n T ) N : Xl < X 2 < ■ ■ ■ < X N }. 



is the section of our hexagon by the vertical line with coordinate t, and 



X NT is the set of all ./V-tuples in this section. 
For any T £ (l(N,T,S) denote by n,T 2 ,. 
numbering starts from the bottom one. Thus, 



tn the corresponding paths, 
1}, 



n = {l - 1 = 7-i(0),7-i(l), • • • ,Ti{T - l),Ti(T) = S + I 

so that n{t + 1) - Ti{t) £ {0, 1}. Note that 

n(t) G %N tT , (n(t), ■ • ■ , Tjv(t)) G <%V iT . 
Consequently, any family of paths T can be identified with a sequence 



{X(1),...,X(T)}, X(t)eX j 



s,t 

AT,T' 



where 

X(t) = (n (£),..., TAr(t)). 

In fact, X(t) is a Markov chain with time t, see [Jl], [J2], |JN| . |Gor| . Its 
transition probabilities are given in Proposition 12 . 21 below. 

Through the Sections [3E] the parameters N and T remain fixed and we omit 

~s,t 

-NT' 



them in different notations. We write 3£ s >* instead of lJ T , X s ' 1 instead of 



X N ' T , and so on. 

In the present paper we introduce a discrete time Markov chain M(r), where 
r is a time variable. M(r) takes values in Q(N, T, r) and one-dimensional dis- 
tributions of M(r) coincide with [i(N, T, r). 

2 Four families of stochastic matrices 

2.1 Properties of sections X(t) 

Denote by ps.t the projection of fJ,(N, T, S) to X s : >*, i.e. 

ps.tij) = Prob{X(t) = Y}, Y = (j/i < • • • < y N ) e X st . 

The following two propositions were proved in |Jli Theorem 4.1] and [Gorl 
Lemma 4]. Below we use the Pochhammer symbol (a)k = T(a + k)/F(a) = 
o(o + l)--(o + A;- 1). 

Proposition 2.1. We have 

N 

psay) = z s ,t n (w - %) 2 n «> s, *(w)> 

l<i<j<Af i=l 

where 

u; s '*(x) 



c!(i + JV - 1 - a:)! (5 + N-1- x)\{T -t-S + x)\ 
and 

_^ (t+l)^ 1 (T-t+l)^ 1 (S-t + N)\(T-S + i-l)\ f t\(T-t )\^ X 
S,t fj (T + l),_ 1 (i-l)l \ T! 

Proposition 2.2. For t = 0, 1, . . . , T - 1, 

Prob{X(t + 1) = (j/i, . . . ,y N ) | X(t) = (»!,..., Zjv)} 

Ilto-Vi)- II (JV + S-Si-l). n fo + T-t-S) 

(t - t) N ■ n fe - ^o 

i/ J/, — Xj G {0, 1} /or every i, and the conditional probability is equal to zero 
otherwise. 



We also need "co-transition probabilities" of the Markov chain X(t). 
Proposition 2.3. For t = 1,2, ... ,T 

Prob{X(t - 1) = ( yi , ...,y N )\ X(t) = (»!,..., x N )} 

XliVj-Vi)- II »i- II (t + N-l-Xi 

i<j i: yi—Xi — 1 i: yi—Xi 



(t) N ■ n (xj - Xi) 

i<j 

if Hi — Xi G { — 1,0} /or every i, and the conditional probability is equal to zero 
otherwise. 

Proof. Straightforward computation using 

Prob{X(t - 1) = ( Vl , . . . ,y N ) | X(t) = (»!,..., sat)} 
= Prob{X(t) = (x 1; . . . , x N ) | X(i - 1) = (y 1; . . . , y N )} ■ Ps^i(yi, ■ ■ ■ ,Vn) 

PS,t( x l> ■ ■ ■ i x N) 

D 

2.2 Stochastic matrices 

We need four families of stochastic matrices P t + , P t i , P s +, PgL- 

Pf^(X,Y) is an \X S ^\ x \X s - t+l \ matrix, X = (xi < ■ ■ ■ < x N ) G X s -\ 
Y=( yi <---<y n )eX s > t+1 ; 

IKvj-Vi) n (N + s-Xi-i) n (»i + r-t-5) 

pS,t/y y\ — iSL i:yi=Xi + l i:yi=Xi 

t+( ' ' (T-tU-Uixi-Xi) 

i<j 

if yi — Xi E {0, 1} for every i, and P 4 _/ (X, V) = otherwise. 

Pg£(X,Y) is an \X S ^\ x |A' S + 1 '*| matrix, X = (zi < • • • < x N ) G X s -\ 
Y=(y 1 <---<y n )eX s+1 < t ; 

Yiivj-vi) n (Jv+t-^-i) n (xi+T-t-s) 

pS,try y\ _ ]<3_ l:yj=Xj + l l: Vi=Xj 

s+[ ' J ~ (T-SU-Uixj-Xi) 

if yi — Xi G {0, 1} for every i, and P S+ {X, Y) = otherwise. 

Pt*{X,Y) is an |;f S ''| x I* 5 '*- 1 ! matrix, X = ( Xl < ■ ■ ■ < x N ) G X s >\ 

Y=(y 1 <---<y n )eX s > t - 1 ; 



ii(yj-yi) n xi n (*+^-i 



Xi 



S,t/\r \ J r\ i <3 i:yi=Xi-l i: yt=Xi 



p t r(x,Y) 



(t) N ■ l\( X 3 ~Xi) 
i<j 



if yi — xt G { — 1, 0} for every i, and P t Z (X, Y) = otherwise. 

P|l*(X,y) is an |A? S '*| x |A" S - M | matrix, X = (x x < ■ ■ ■ < X N ) € X S '\ 
Y=(y 1 <---<y n )eX s - 1 > t ; 

n(w-w) n ^ n (s+jv-i-so 

pS,t(y y\ _ ]<0_ i:yi=Xi-l l:yi=Xj 

S - ( ' '~ (S)N-Il(Xj-Xi) 

i<j 

if yi — xt G { — 1, 0} for every i, and P s '_ (X, Y) = otherwise. 

Looking at spaces that parameterize rows and columns of these matrices one 
can say that P t + increases t, P t l decreases t, while P s ^_ increases S and P s l_ 
decreases 5. 

Theorem 2.4. All four types of matrices defined above are stochastic, and they 
preserve the family of measures ps,t- I n other words 

]T P t ¥(X,Y) = l, J2 Pt±{X,Y) = l, (1) 



Ps,t±i(Y)= £ P t ±(X,Y>PS,t(X),PS±i,t(Y)= E Ps±{X,Y)-p s . t {X). 
xex s - f xex s - f 

Proof. Propositions HH and E3] imply the claim for P£f(X, Y) and Pfl^X, Y). 
Now observe that the space X '* is unaffected when we interchange param- 
eters t and S, i.e. 

Moreover, the measures ps,t are also invariant under S «-*• t, i.e. 

Psj = Pt,s- 

Finally, note that Pf*(X,Y) becomes Pg£(X,Y) under S <-» t and 
P^X.F) becomes Pg'*(X,Y). 

Consequently, applying S «-» t to the relations for P t ± we obtain needed 
relations for P$± . □ 

2.3 Determinantal representation 

In this section we write our stochastic matrices in a determinantal form. This 
representation is very convenient for various computations. 
First, we introduce 4 new two-diagonal matrices: 

{N + S-l-x, i£y = x + l, 
T-t-S + x, tiy = x, a; G X s '*, y£X s > t+1 ; 

0, otherwise, 




x e X 5 '*, y e X 5+M ; 



otherwise, 



X, if y = x — 1, 

Utl\x,y) = { t + N-l-x, iiy = x, x e X s '*, y e X s '*" 1 ; 

0, otherwise, 



if y = x-l, 
Ugt{x,y) = { S + N-l-x, iiy = x, x e X s '*, y e X s " 1 '*. 

0, otherwise, 

It is possible to express stochastic matrices P t ± , P§± through certain minors 
of the matrices defined above. 



Proposition 2.5. We have 



n (vj - Vi) dct P t + fa, yj)]i,j=i,...,N 



p£«(X,y) = !5L 



[t - *)n • n (»j - *») 

11 (Vj - Vi) det[f/ s ^(xi, j/j)]i,j=i,...,jv 



Ps + (x ' y) - ' (T-^-nfe-^) 






II (W - y») det[[/ t J. (xi,3/j)]i,j=i,...,iv 



p t s :*(x,y) = ^- 



(*)jv • II ( x i - ar») 
II (% - I/i)det[C/ s ^(a;i,j/j)]i 1 j=i,... 1 jv 



p|i*(x,y) = — 



(s)n ■ n ( x j - x i 

i<j 



Proof. Straightforward computation using the definitions of stochastic matrices 
P t ± , Pg± and matrices U t ± , U s ±. 

Any submatrix of a two-diagonal matrix, which has a nonzero determinant, 
is block-diagonal, where each block is either upper or lower triangular matrix. 
Thus, any nonzero minor is a product of suitable matrix elements. □ 



2.4 Spectral decomposition of stochastic matrices 

In this section we modify the determinantal representation of the stochastic 
matrices and introduce new representation, which we call spectral decomposi- 
tion. Spectral decomposition is of crucial importance for computing correlation 
functions of the processes that will be constructed later on, and for finding bulk 
limits of the processes. Results of this section will be used in Sections \5\ and [6l 
while Sections [3] and Q] are independent of these results. 

Let us introduce some notations. 

Denote by H k ' (x) the Hahn polynomial of the degree k corresponding to the 
parameters N, T, t, S. Domain of definition of these polynomials coincides with 
X '*, and the polynomials are orthogonal with respect to the weight function 
w '*(x) defined in Proposition ^. 11 For definition and explicit formulas for Hahn 
polynomials see |KSj . more information about the usage of Hahn polynomials 
in our case can be found in [Gorj . 

Let us denote by iff k ' (x) the normalized Hahn polynomials 

Here (H k ' , H k ' ) is the squared norm of the polynomial H k ' in 
L2 (3~ S '', w s,t (x)). The functions ^ k ' (x) form an orthonormal basis in the space 
L2(X s,t ) (this Li is with respect to the uniform measure on 3t '*). 
Let 



t + N J V T + N-t- \ 

Note that c t _J_ (i) does not actually depend on S, but it is convenient to use it 
for the notation. 

Finally, denote by v t ? the |X S '*| x |X t+1,s | matrix given by 

i>0 

The following proposition was proved in Gor, proof of Proposition 5] 

Proposition 2.6. Let X — {x\ < ■ ■ ■ < Xn) and Y — (j/j < ■ ■ • < y^) be 
elements of X '* and X ' t+1 , respectively. Then 

pS,tf Y v\ - V / Ps,t+i(Y)det[vf^(x i ,y j )] t ^ h .^ N 
^t+ V A , Y ) - — 1 , j^i 

v^w n 4*«) 

Stochastic matrices P t J_ admit similar spectral decomposition with trans- 
posed matrices. 
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Denote 

s,t s,t-i 
c t- -=Ct+ 



and 



S.t / S,t-1\T 

v t- : = {V t + ) 



(here ( ) T means matrix transposition). 

Proposition 2.7. Let X = (an < ■ ■ ■ < xn) and Y — (yi < ■ ■ ■ < yjj) be 
elements of X ^ and X s,t ~ l , respectively. Then 



T}S,t fY v s _ Vps,t-i(Y) det[v t l {x u yj)) iii=1 ,... tN 
P t _(X,Y) = - == ^-— 

n cfi*(») 



i=0 



Proof. Recall that while P t ■! coincides with transition matrix of the Markov 
process X(t), matrix elements of P t l are co-transition probabilities of Propo- 



sition [273] Consequently, 



?S,t lv v-v _ AS.t-lQO r>S,t-l 

psAx) 



p?i\x,y) = P5it ~ l1 j • p^-^yx) 



/t>s,t_i(F) y/p s ,t( x ) det[«4* : (yi)»j)]i,j=i,.... 



A r 



<*,(*) y^T-TIF) ^^W 



( = 

5,t/ 



\//35,t-i(y) det[u t I (yi,Xj))jj =1 ,... !N 

i=0 



D 



Define 



4i(0 = i/(l-o-htr)(l- 



" s+w yv s+n/ \ t + n-s-i 

(these constants do not acually depend on t). Set 

i>0 

Recall that matrices P t _£ and P s l are connected by the involution t ■*-» S, 
thus, Propositions 12.61 and 12. 71 imply the following statements. 

Proposition 2.8. Let X = (an < • • • < sat) emd Y" = (j/i < • • • < J/jv) ^e 
elements of X ^ and X s+lt , respectively. Then 

p3,t fY V\ - y/PS+lA Y ) dct [ v s'U x i>y3)]i,3=l,-,N 

v^w n 4i(i) 

»=0 
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Proposition 2.9. Let X = (x\ < ■ ■ ■ < Xn) and Y = (y\ < ■ ■ ■ < yjv) be 

elements of X s,t and X s ~ l ' f , respectively. Then 

p s.t (x y) = VP5-i,t(y)det[z;f: 1 ' t (y i ,x i )] i , j ^i„„,Af 



n 4<ni) 

i=0 



Theorem 2.10. The families of stochastic matrices P t ± and P s ± commute, 



2.5 Commutativity 

Theon 

that is 



r>S,t pS,t+l _ r,S,t pS-l,t 
r t+ ' "s- ~ r S- ' ~t+ ' 

**- ' "S- ~ ^S- ' **- ' 

pS,t pS,t+l _ pSi p S+l,f 
**+ ' ^5+ — r S+ ' / t+ ' 

pSi pS,t-l _ p Si pS+l.f 
**- ' ~S+ ~ *S+ ' "t- ' 

/or any meaningful values of S and t. 

Proof. Proofs of all four cases are very similar and we consider only the first 
one. 

{p s,t . P |,*+i )(X) y) = Y, P t S +( X > Z ) ■ P s- + \Z, Y) 
ze*s,t+i 

I\(Vi -Vi) E det[t/ tH l*(a;i,2;j)]i 1 j=i,...,jvdet[t/' s ^ +1 (^,3/j)]i,j=i ) ...,iv 
(T - t)jv • (S')jv • EI (»i ~ »j) 

Applying Cauchy-Binct identity we obtain 

^ det[!7 tH l*(x,, z,-)]i,j=i,...,jv det[C/ s l* +1 (zj, %-)]j,j=i,...,jv 
zeArs,t+i 

= det [([/(/ • C/ s l* +1 )(a;i,yj)]i,j=i,...,Jv- 
Thus, 

r S,t r rS,t+l\ 



(P^*.p|^+ 1 )(jf,y) = ^- 



II (Vi - Vi) det i( u t+ ■ U s '_ )(xi,yj)]i tj =i t ... t N 



(t - t) N ■ (s) N ■ n ( x i - x j) 

i>3 



Similarly 



jS.t nS-l,tw V 1^\ »>J 



(j/i - yj) dct i(Us- ■ U? + 1 '*)(a;i, 3/j)]i,j=i,--,Jv 



(p^-p^XXT) 



(T - t)jv • (5);v • II ( x i _ x j) 

i>3 
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Our problem reduces to verifying the equality 



where 



T rS,t TT S,t+l rrS,t TT S-l,t 
U t+ ' U S- ~ U S- ' U t+ 




,rd computation one proves that 




u t+s- — u t+ • u s- ~ u s- ' u t+ ' 




f (N + S - 1 - x) (N + S - 2 - x), 
\ (N + S - 1 - x)(T - t - S + 2x + 1), 
\x(T-t-S + x), 
0, 


if y — x + 1, 
if y = x, 
if y = x- 1, 
otherwise. 



E#* s _ (*>») = 



D 

Remark. Another way to prove the commutativity is to use the spectral 
decomposition introduced in Section [2.41 and to observe that coefficients c t _j_ (i) 
do not depend on S, while coefficients c s '±(i) do not depend on t. One computes 



w s,t {x)/w s-x,t + i {y) 



(t + N)(T + N -t - 1)(S + N - 1)(T + N - S) 

Es,tf\ s,t+i/.\, T .s,t/ \ lTr s-i,t+i/ -, 



ES,tf\ S-l,t/.\. T ,S,t/ \, T ,S-l,t+l/ \ 



i>0 

3 Markov step 5 i-> S ± 1 

In this section we aim to define two new stochastic matrices 

P# MS+1 (X,F), X(=n(N,T,S), Yen(N,T,S + l) 

and 

^fLs-iC^F), XeQ(N,T,S), Y €ti(N,T,S-l) 

that preserve the measures [/,(N,T,S). Both -PfLs+i an d P§i-+S— l depend on 
parameters AT and T but we again omit these indices. In Introduction we called 
these matrices -Pf xfcxc with a = T — S, b = S, c = N. 

Suppose we are given a sequence X= (X(0), X(l), . . . , X(T)) G n{N,T,S) 
(recall that X(t) G A" 5 '*). Below we construct a random sequence Y = 
(Y(0), . . . , Y(T)) G tt(N, T, S + 1) and therefore define the transition proba- 
bility (or, equivalently, stochastic matrix) Pg^s +1 (X, Y). 

First note that F(0) G A' s+1 ' and \X S+1 >°\ = 1. Thus, F(0) is uniquely 
defined. We will perform a sequential update. Suppose Y(0), Y(l), . . . , Y(t) 
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have been already denned. Define conditional distribution of Y(t + 1) given X, 
Y(0),F(l),...,Y(t)by 



Prob{F(i+l) = Z} = 



(2) 



p t + +1,t (y(i), z) • pf + M+1 (z, x(t + 1 )) 

i pS+i-t p s+i,t + i ){Y{thx{t + 1)) 

_ P^ +1 (X{t + 1),Z)- P t s+M+1 (Z, Y(t)) 
(P s ^ +1 Pl +1 - t+1 )(X(t + l),Y(t)) 

(The second equality follows from ps+i,t+i{X)P t ^~ ' t+ (X,Y) = 
Ps+ht (Y)P t S + +1 - t (Y,X).) 

This definitions follows the idea of [DF[ Section 2.3], see also [BF . 

Observe that (p£ +M Pf+ M+1 )(y(f),X(£ + 1)) > 0. Indeed 

(i& +1 -*p#± M+1 )(y(t),x(t + i)) = (p^ 1 - t p t s + t )(Y(t),x(t + 1)) 

> P$W(Y(t),X{t)) •P t S /(I(t) ) I(t + l)) > 0, (3) 

because F(i) was chosen on the previous step so that P s ^ ' (Y(t),X(t)) > 0. 

One could say that we choose Y(t+ 1) using conditional distribution of the 
middle point in the successive application of P t+ ' and P s _ ' (or Pg\ 

and P t _ ), provided that we start at Y(t) and finish at X(t + 1) (or start 

at X(t + 1) and finish at Y(t)). 

After performing T updates we obtain the sequence Y. 

Equivalent^, define PjL s+1 by (cf. [EH Section 2.2]) 

"r-i p£ +1 '*(y(t), Y(t + 1)) • pf + M+1 (r(t + i),x(t + 1)) 



P^ S+1 (X,Y) 



' ' /T,S+1,« d s+m+i 



(i^^P|+^^)(y(t),X(t + l)) 



t+ ^ s- 
if T fl\p t s t 1 ' t Pl± 1 ' t+1 )(Y(t),X(t+ 1)) > 0, 

0, otherwise. 



Theorem 3.1. TTie matrix PJL S+1 on n(JV, T, S) x ft (AT, T, S + 1) is stochas- 
tic. The transition probabilities Pg^g, 1 (X, Y) preserve the uniform measures 
H(N,T,S): 

fx(N,T,S + l)(Y)= J2 P^ s+1 (X,Y)n(N,T,S)(X). 

xen(N,T,s) 

Proof. First, let us prove that the matrix Pg^s+i ^ s stochastic, equivalently: 
Yl P t T-\Y(t),Y{t + 1)) ■ P s s + l - t+ \Y(t + l),X(t + 1)) 

Y(0),...,Y(T) t=0 ( P t+ ' P S- ' )( Y ( t )> X ( t + !)) 
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where the summation goes over all (Y(0), F(l), . . . , Y(T)) e ti(N, T, 5+1) such 
that 

T-l 

J] (pf + +1 - t p|+ M+1 )(r(t), x(* + 1)) > o. (5) 

t=o 
We write (j4]) in the form 



e n 



F f ; +1 '(y(t), y(* + 1)) ■ p|+ • +i (r(i + 1), x(t + 1)) 



Y(0),...,Y(T-1) t=0 ( P *+ ' P S- )( Y ( t )' X (t + l)) 



X 



p 5+i,T-i (y(T _ 1); y(T)) . P f +i.r (y(T)) X(T)) 

zL, (p i s : +1 ' T " 1 p|+ 1 ' T )(r(T-i),x(r)) 



y(T) *> *+ 

Summing over Y(T) we obtain 

y. T-3 f S + i, t(y(f)| y(f + 1)} . p|+i.*+i (y(f + i^yft + x)) 

y(o),...,y(T-2) *=o (-^4- ' p s- ' )(^(*)>^(*+ 1)) 

p g+i,T-2^ y ^ T _ 2),y(T - 1)) ■ P^^-^P - 1), JC(T - 1)) 

X yfcLq ( p5+l,T- 2p |+l,T-l )(r(r _ 2)) X{T _ 1)} ■ 

Next, we want to sum over Y(P — 1). Inequality ([5]) implies that the summation 
goes over Y(T - 1) such that (P t + +1 ' T ~ 1 P|+ 1 ' T )(r(T - 1),X(T)) > 0. Note 
that 

(p 5+l,T-l p J+l,T )(r(T _ 1))X(T)) = {p S+l,T-l p S,T-l )( Y {T lhX{T)) 

> P^ 1 ' T - 1 (Y(T-l),X(T-l))P t s f- 1 (X(t-l),X(T)). 

Consequently, if (P t + +1 ' T_1 Pf± 1,T )(Y(T - 1), X(T)) vanishes, then 
P s _ ' ~ (Y(T — 1), X(T — 1)) vanishes too. Thus, we may drop the inequality 
that restricts the summation and sum over all possible Y(T — 1). After that we 
sum over Y(T — 2), and so on. After summing over all Y(t), t — T, T — 1, . . . , 1 
and noticing that Y(0) has just one possible value we arrive at (TJ|. 

Now we prove that the transition probabilities Pg^g +1 (X, Y) preserve the 
uniform measures /j,(N,T, S). It is equivalent to 

»(N,T,S+1)(Y)= £ pLs+i(X,YMN,T,S)(X). (6) 

X = (X(0),X(1),...,X(T)) 

Since X(t) can be viewed as a Markov chain with time t, 

v(N,T,S)(X) = m° s (X(0))-P t s > o (X(0),X(l))...P t s l T - 1 (X(T-l),X(T)), 
where m l g(X(0)) is the unique probability measure on singleton X S: °. 
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Thus, the right-hand side of © is equal to 



T-l 



]T m° s (X(0)) J] if/(X(i),X(t + l)) 

X(0),...,X(T) 



t=0 

T-i pg+i.«(y (f)i y (f + 1)) . P S+^+\ Y (t + 1), X(t + 1)) 



t=o 



(p 5 + l, tp S+l, t+ l )(y(i)5X(t + 1)) 

Note that 

m°(X(0)) = m° s+1 (r(0)) = i^+^CFfO), Jf(0)) - P^ T (Y(T),X(T)) = 1 
and write ([7]) in the form 

T-l 



m 



s+ 



1 (Y(0))]JPf+ 1 ' t (Y(t),Y(t + l)) 



t=o 



x Z^ 11 /d s+m d s,* n 



X(0),...,X(T) t=0 



(P$™Pg)(Y(t),X(t + l)) 



(We used the equality Pfj" 1 '^' = p£ +1 '*Pfj' M+1 .) Summing first over X(0), 
then over X(l), and so on, we get 



T-l 



m° s+1 (Y(0)) J] P^ +1 '*CK(*),y(t+l)), 

which is exactly the distribution n(N, T, S + 1)(Y). 

Similarly to Pg^s+i) onc defines a transition matrix 

p|^ s _ 1 (x,r), xen(N,r,s), y gQ(n,t,s-i), 



n 



by 



' T — 1 T>S— l»t 



pJL s -i(x,y)H 



n 



P.V 1 ' 1 to, n* + 1)) ■ p#+ M +1 (n* + i), *(* + 1)) 
(p&- M p#+ x ' t+1 )TO,*(* + i)) 



T-l 



t=0 

0, otherwise. 



Similarly to © there is another way to write P§ l _^ s _ 1 because of the equality 



■>s-i,t 



p t V^(Y(t),y(t + i))-p 



S-l,t+l 



(Y(t + l),X(t+l)) 



{ pS-l,tpS-l,t + l )( y { ^ X{t+1)) 

_ P^ +1 (X(t + 1), Y(t + 1)) ■ P t s _-^ t+1 (Y(t + 1), Y(tj) 
(P|l* +1 P t i~ 1,t+1 )pr(i + 1), K(t)) 
Similarly to Theorem 13.11 one proves the following claim. 
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Theorem 3.2. The matrix P^s-i on fi W T,S) xfl(N,T,S -I) is stochas- 
tic. The transition probabilities Pg l _ ¥ g_i(X,Y) preserve the uniform measures 

V(N,T,S): 

H(N,T,S-1)(Y)= J2 PLs-i(^Y)KN,T,S)(X). 

XgQ(N,T,S) 

Remark. The above construction performs sequential update from t = to 
t = T . One can equally well update from t = T to t — by suitably modifying 
the definitions. The resulting Markov chains also preserve the uniform measure 
/i(iV, T, S), and they are different from the Markov chains defined above. 

4 Algorithmic description 

In this section we suggest an algorithmic description of the Markov chain from 
the previous section. 

Denote by D(a, b, n) the probability distribution on {0, 1, . . . , n} given by 



Prob({/c}) = D(a,b,n){k} = 



{a)k 



j=0,...,n 



_ a(a + l)...(a + k-l)(b + k)...(b + n-l) 

~£;u a («+ i )-->+.7- i )( & +.7)---(&+«-ir 

4.1 Algorithm for S i— > S + 1 step. 

Suppose we are given X = (X(0),X(1),...,X(T)) E Q(N,T,S). We want to 
construct Y = (Y(0), F(l), . . . , T(T)) 6 0(iV, T, 5 + I). 
In the first place we note that Y(0) is uniquely defined, 

F(0) = (0,l,...,iV-l). 

Then we perform T sequential updates, i.e. for t = 0, 1, . . . T — 1 we construct 
Y(t + 1) using Y(t) and X(t + 1). Let us describe each step. 

Let Y(t) = (j/i < y 2 < •• • < j/jv) and X(t + 1) = («i < x 2 < ■ ■ • < *jv)- We 
are going to construct Y(t + 1) = (zi < z 2 < ■ ■ ■ < z N ). 

Recall that 

Zi G X s+1 ' t+1 = {x e Z | max(0, t + S-T + 2)<x< min(t + TV, 5 + 2V)}. 

T(t) and X(t + I) satisfy ([3]). This implies that Xi — yi is equal to either — 1, 
or 1 for every i. 

• First, consider all indices i such that Xi — yi = 1. For every such i we set 

Zi Xi . 

• Second, consider all indices i such that Xi — yi = — 1 and set Zj = j/j. 
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• Finally, consider all remaining indices, i.e. all i such that Xi = j/». Divide 
the corresponding x^s into blocks of neighboring integers of distance at least 
one from each other. Call such a block a (k, Z)-block, where k is the smallest 
number in the block and I is its size. Thus, we have 



x l =y l = k, x i+1 = y i+1 = k + 1, 



Xt+i-i = Vi+i-x = k + I - 1 



and 



3/*— 1 <k-l, y l+ i > k + l. 

For each (fc, Z)-block we perform the following procedure: consider ran- 
dom variable £ distributed according to D(k + T — t — S — l,k + 1,1) (£'s cor- 
responding to different (fc,Z)-blocks are independent). Set Zi = Xi for the first 
£ integers of the block (their coordinates are fc, fc + l,...,k + £ — 1) and set 
Zi = Xi + 1 for the rest of the block. 

At Figure 4 we provide an example of constructing Y(t + 1) using X(t+ 1) 
and Y(t): there is only one (k, Z)-block and it splits into two groups, here £ = 2. 



X(t) X(t+1) 



Y(t) Y(t+1) 



M 
u 

I 




-*- could not jump (first case) 



^_ split point detemiined by (8) 
(block case) 



*- was forced to jump (second case) 



S S+l 

Figure 4. Example of (fc, I) -block split, I = 4, £ = 2. 

Theorem 4.1. T/ie algorithm described above is precisely S ^ S + 1 Markov 
step given by -PjL+s+i- 

Proof. As was shown in the previous section, the transition S ^ 5 + 1 consists 
of T updates. Namely, given Y(t) and X(i + 1) we define Y(t + 1) by 

P?, +1>t (Y(t). Z\ ■ p s + 1 ' t + 1 > 
Prob{F(t + l) = Z} 



i + 



\Y{t),Z).P s s + 1 - t+l (Z,X{t+l)) 



(p S ++ M jP |+M+i )(rWjX(t + 1)) 

Let Y(t) — (yi < y 2 < ■ ■ ■ < J/at), X(t + 1) = (xi < x 2 < • • • < x N ) and we are 
defining Y(t + 1) = (z\ < z^ < ■ ■ ■ < Zn). Inequality ([3]) implies that for every 
i the difference Xi — yi is equal to either — 1, or 1. Thus, we have 3 cases: 
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1. Xi — yi = 1: Definitions of P t _J ' and P s _ ' imply that both z% — yi 
and Zi — Xi must be equal to either or 1. Consequently, z% is uniquely 
defined and z% = Xi- 

2. Xi — yi = — 1: Again % is uniquely defined, and z% = j/». 

3. Xi — j/i = 0: This is the only nontrivial case. Here one has two possibilities, 
either z% = Xi or z% = Xi + 1. 

When we pass from Xj to Zi, every {k, I) -block is split into at most 
two groups. Namely, we have I + 1 possibilities for the split point 
j € {0,1, ... ,1} 0, 1 . . . , I: Zi — Xi for the lowest j points of the block 
and z-i = Xi + 1 for the other points of this block. 

Now we want to compute the probabilities of different splits of the blocks. 
We have 



Prob{y(t+l) = Z} = 



P t S + +ht (Y(t), Z) ■ Pf + M+1 (Z, X(t + 1)) 
(p s +Mp s+ M +i )(mx(t + 1)) 



= n (vi+T-t-s-i)- n w+s-vi) 

i-.Zi—y^ i:Zi—yi-\-l 

x TT (N + S - Xi) ■ TT (xi + 1) • (factors independent of Z) 

i:Zi—Xi i:Zi—Xi-\-l 

This formula implies that blocks split independently For each (k, Z)-block 
the probability of split position j is equal to 

k+j-i k+l-i 

]J ( a + T -t- S -1)(N + S -a) • JJ (a+l){N + S-a) 

a—k a—k+j 

x (factors independent of j) 

Since (A^ + S — a) is present in both products, this probability can be 
written as 

k+j-l k+l-l 

Y[ (a + T-t-S-l)- Yl (a + 1) • (factors independent of j) 

a—k a—k+j 

which is exactly the distribution D(k + T — t — S — 1, k + 1, 1). 

□ 

4.2 Algorithm for S i— ► S — 1 step 

The S <-^ S — 1 step algorithm is very similar to the S ^ S + 1 one. 

Suppose we are given X = (X(0),X(1), . . .,X(T)) e Cl(N,T,S). We want 
to construct Y = (Y(0), Y(l), . . . , Y(T)) e Sl(N, T, S - 1). 
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As above, note that Y(Q) is uniquely defined, 

y(o) = (o,i,...,Ar-i). 

Then we again perform T sequential updates, i.e. for t = 0, 1, . . .T — 1 we 
construct Y(t + 1) using Y(t) and X(t + 1). Let us describe each step. 

Let Y(t) = (yi < y 2 < ■ ■ ■ < Un) and X(t + 1) = («i < X 2 < ■ • ■ < x N ). We 
are going to construct Y(t + 1) = (z\ < z 2 < • ■ • < Zn)- 

Recall that 

Zi e X 5 - M+1 = {x e Z | max(0, i + S - T) < x < min(i + N, S + N - 2)}. 

Y(t) and X(t+1) satisfy (P t + _1, *P| 4 : 1 '* +1 )(y(t),X(tH-l)) > 0. This implies 
that Xi — yi is equal to either 0, 1 or 2 for every i. 

• First, consider all indices i such that Xi — yi — 0. For every such i we set 

Z% Xi . 

• Second, consider all indices i such that Xi — y, = 2 and set Zj = j/j + 1. 

• Finally, consider all remaining indices, i.e. all i such that Xi = yi + 1. 
Divide the corresponding x^'s into blocks of neighboring integers of distance 
at least one from each other. Call such a block a (k, l)'-b\ock, where k is the 
smallest number in the block and I is its size. Thus, we have 

x l = y l + l = k, x i+ x = j/j+i + 1 = k + 1, ..., Xi+i-i=yi+i-i=k + l — l. 

For each (fc, Z)'-block we perform the following procedure: consider random 
variable £ distributed according to D(N + t — k + l,N + S — k — 1,1) (^'s cor- 
responding to different (k, //-blocks are independent). Set Zi — yi for the first 
^ integers of the block (their coordinates are fc — l,k, . . . ,k + £ — 2) and set 
z i = Vi + 1 for the rest of the block. 

Theorem 4.2. T/ie algorithm described above is precisely S i— ► S — 1 Markov 
step defined by Pg,_ > c_ 1 . 

The proof is similar to Theorem 14. II and we omit it. 

4.3 Numeric experiments 

The S i— ► S ± 1 steps can be used in different ways. 

1. Suppose that our aim is to sample a random tiling (equivalently, random 
family of paths) T € H(7V, T, S) from the uniform measure fi(N, T, S). 

We start from the unique family of paths % 6 0(JV, T, 0). Indeed, 
|fi(iV,T,0)| = l. 

Next, we perform 5 steps. During rth step we construct 7" r G fl(N,T,r) 
distributed as fi(N,T,r), using already constructed family of paths T r ^\ £ 
£l(N,T,r — 1). Theorem 13.11 implies that Tg is the desired random element 
of n(N,T,S). 

Let us estimate the number of operations. Every update takes O(N) oper- 
ations. For every S — > S + 1 step we have to perform T updates. Consequently 
to sample from il(N,T,S) we need O(NTS) operations. 
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On Figure 5 we show a random surface generated by our algorithm. Here 
iV = 1000, T = 2000, S = 1000. It took less than 4 minutes on our laptop (Intel 
Core2 Duo 2.2GHz, 2Gb Ram) to generate this tiling. Theoretically predicted 
"arctic ellipse", see |CLPj . is clearly visible on our picture. 





Figure 5. Random surface corresponding to the big tiling. 

2. Using our steps one can also construct equilibrium dynamics S i— » S+l 
S or S h-» S — 1 h- ► S. 

On Figure 6 we show the evolution of the "filled box" tiling under S 
S+l^ S dynamics. Here N = 50, T = 50, S = 20. 
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Original tiling 



After 20 steps 




After 100 steps 




After 1000 steps 




Figure 6. Evolution of tiling under S <—* S + 1 i— ► S dynamics. 

5 General 2-dimensional dynamics, its sections 
and correlation functions 

5.1 Construction of dynamics and its correlation functions 

In this section we construct a family of Markov chains on the spaces Q(N, T, S) 
using the transition probabilities Pg^g +1 (X, Y) and Pg l ^s_ 1 (X, Y) introduced 
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in Section 3. Namely, we want to combine 5^5+1 and S i— ► S — 1 steps. To 
fix the order of "+1" and "—1" steps we introduce an auxiliary sequence {e;} 
of +l's and -l's. 

Formally, let < So < T and let e = {ej}i=i,2,... be an arbitrary finite or 
infinite sequence of + l's and —l's such that 



0<S + J2 e ^ T 



n=l 

for every m. The last condition is necessary to ensure that the state spaces 
Q(N, T, S) of the process are not empty. 
For any integer r we denote 

Given So and e* let us define a Markov chain 
9K So ,e(r), r = 0,l,... 
3^S ,e( r ) takes values in fl(N,T,S(r)), its initial distribution is fi(N,T, So)- 

Prob{Af Soi£ (0) = X} = M (iV,T,Sb)({X}). 
Transition probabilities of our process are given by 

Proh{m So Ar + l)=Y I 9%,, e (r) = X} = (Sf +1 ^' ^' * £ " +1 = * , 

Theorems 13.11 and 13.21 imply that one-time distributions of 9Jts 0j£ (r) are 
exactly fi(N,T,S(r)). 

Example 1. If So = and {e^} = {1,1,1,...}, then 9Jts ,e(Y) is precisely 
the chain used for the random tiling sampling in Section [4.31 

Example 2. If we set e, = (—1)* and then restrict dJls ,e(r) on even r we 
get a stationary Markov chain from Section 14.31 

Recall that X G Q(N,T,S) is a sequence X = (X (0),X(1), .. . ,X(T)). 
Given a trajectory of the Markov chain 9Jts , e (r) we can construct a point 
configuration in 1? with coordinates (r, t, x) such that the point (ro,to,xo) is 
occupied if and only if xq £ (9Jts Oje (r )) (t ). Thus, our Markov chain defines 
a measure on such point configuration or, equivalently, a random point process 
in Z 3 . Denote it by M. 

Define the nth correlation function of M by 

R n (ri,ti,x 1 ;r 2 ,t2,x 2 ; ■ ■ ■ ;r n ,t n ,x n ) 
= M{M e Conf(Z 3 ) | (n,ti,xi) eM,(r 2 ,t 2 ,x 2 ) G M, .. . , (r n ,t n ,x n ) £ M}. 

These correlation functions uniquely define the process M. 
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Through the rest of the paper we concentrate on computation of correlation 
functions R n . Unfortunately, we are not able to fully describe R„ for all possible 
arguments. But we can compute R n on certain two-dimensional sections of Z 3 . 

The main result of this section is the following statement. 

Theorem 5.1. Let r\ < r 2 < ■ ■ ■ < r n , t\ > ti > • • • > t n . Then 

R n (r 1 ,t 1 ,xi;r 2 ,t2,x 2 ;- . . ; r„, t n , x n ) = det [K{r i ,t i ,Xi;rj,tj,Xj)] i ^ =x >n , 



where 



N-l 

K(r,t,x;r',t',x') = £ _^*fM.* (a; )*fM.* ( tf ), 

if r > r ,t <t i 
K(r, t, x; r',t',x') = -J2 ^' t '^ r) '\x)^ ) '^ (y), 

i>N 

if r < r or r = r ,t > t ; 



C/ = 1, 



r'-l 



<A v, * = -q c s(r )<k{i) -q c 5(,), t(f)) 

fc=t' + l fc=r 

where Sek+i stands for S+ if €k+i = +1 and S— otherwise. 

In Section[5]we will study the bulk asymptotics of these correlation functions. 
For n = T2 = • • • = r n , Theorem 15.11 was obtained in [J2], |JN] . [Gorj . 

5.2 Admissible sections 

We call a sequence A = ((ro,to), (ri,ti), . . . (r n , £„)) an admissible section of Z 
provided that 

1. (r i ,i i )e{0,l...}x{0,l,...,T} 

2. r < ri < • • • < r n 

3. to > tx > ■ • • > t„ 

4. For every i = 0, 1, . . . ,n — 1 either r-j+i = r*j + 1, £,;+i = ti or r^+i = ri, 
ti+i = U — 1 

Figure 6 gives an example of an admissible section. 
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Figure 6. Admissible section. 

Given an admissible section A we introduce a sectional process M^/i), h ■ 
0,1..., n, by 

M A (h) = (Wls , € (rh))(th) 

Theorem 5.2. M. A (h) is a Markov chain with initial distribution 

Prob{M A (0) = X} = p S N ^ )Ao (X) 
and transition probabilities 



Prob{M A (h + 1) = Y | M A (h) = X} 

■pS(r h ),t h{xX) ^ ifth+i 



th - 1, 



= <P 



S(r h ),t h 
S+ 
yS(r h ),t h 



(X,Y), ifS(r h+1 ) = S(r h ) + l, 
(X,Y), ifS(r h+1 ) = S(r h )-l. 



Proof. This theorem follows from [BF| Proposition 2.7]. Let us explain the cor- 
respondence between our notations and notations of |BFj . To avoid confusions 
we denote by r the time variable t from |BFj . Then r corresponds to time 
r of Markov chain 97ts , e (r); k of |BF] corresponds to t; state space <Sfc(V) is 
X S(r),k. matrices Aj^_ 1 (- , • | t) of HE) are p^' k ( . , . ) ; matrices P fc ( ■ , ■ | r) 
are PgJ (•, -)i the commutation relations Aj^P/j-i = -PfcA^j _ 1 are exactly 
PgfPf.^ 6 ' 1 — P t _lPgl~ above; t is r and m„(-) of |BF] corresponds to 

PS(r ),t - n 

5.3 Correlation functions 

To compute the correlation functions R n we are going to use a variant of the 
Eynard-Metha theorem (see |EMj and |BO| Section 7.4]). Let us state it first. 
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Proposition 5.3. Assume that for every time moment h we are given an or- 
thonormal system {g^} n >o ^ n linear space ^({0, 1, . . . , L}) and a set of numbers 
Cg, cj, . . . . Denote 



\,h+l{x,y) = J2 C n9n(*)9n +1 (y)- 



I'h 



Assume also that we are given a discrete time Markov process Vh taking 
values in N -tuples of elements of the set {0, 1, . . . , L}, with one-dimensional 
distributions 



, 2 
N, 



(det[^_i(si)] w=li 

and transition probabilities 

det [vh,h+i(xi,Vj)] iiJ=lt ... tN teb [9i-i(Vj)] iij=1} ... tN 



det[<#_i(*i)] ii=1 N U 4 
Then 

Probjxi e P kl ,...,x n e Pfc n } = det [if(fc i ,x l ;fc :/ ,a; i )] 4j=:L 



wher 



N-l 



K{k,x;l,y) = Y, —gH x )9l(y), k>l; 

i=0 C i 

K(k, X] l,y) = -J2 c^gH^gliy), k < I; 

i>N 

^ = 1,^=^.^ ci" 1 . 

Now set P/i := M^(/i). Then we can take orthonormal functions g\ to be 
the functions ^ (rh) '* h (iz;) defined in Section H3 and 

(ct rh) ' th {^}, if *fc+l = t& - 1, 

c? := ^ cf£ h) ' th (*), if ^(rfc+i) = S(r h ) + 1, (9) 

\cf: h) ' th (i), HS(r h+l ) = S(r h )-l, 

(vfl rh) ' th (x,y), i£t h+1 =t h -l, 
v Kh+1 (x,y) := } v S s ( l h) ' th {x,y), if 5(r h +i) = S{r h ) + 1, 
Uf- M, ifS(rfc+i)=5(r fc )-l. 

Proposition 5.4. Markov chain M^(/i) satisfies the assumptions of Proposi- 
tion\5.3[ 



26 



Proof. By Theorem 12.41 one-dimensional distributions of JA A (h) are given by 
the measures Ps(r h ),t h - These are Hahn orthogonal polynomial ensembles (see 
e.g. [j2]. [Gorj ). It is well known (see for example [Konl Section 2.7]) that this 
distribution can be written in the form 



Ps(r h) ,t h (X) = (det[^ r 1 ^(x j ) 



2 
i,j=l,...,N, 



Propositions 12 . 7H2 .91 imply that the transition probabilities can be expressed in 
the required form too. □ 



Applying Proposition 15.31 we obtain the following 
Proposition 5.5. 

Probjiri £ M A {hi),.. . ,»„ G M A (h n )} = det [K (k t , x t ; kj , x 3 )] ij=1 ^ n . 



N-l 



K(k,x;l,y) = £ 4z*?™' th ( x )*?™' tl (.v)> k * '? 

z=0 C i 

K{k, x;l,y) = -^ c*'*f ^ (z)*f W '« {y),k< I; 

i>N 
k,k _ 1 k,l _ k fc+1 J-l 

and coefficients c[ are given by (0). 

Proof of Theorem \5.1l If ri < r2 < • • • < r n , t\ > ti > • • • > £„, then the 
sequence {(ri, ti), . . . , (r n ,t n )} can be included into some admissible section A. 
Applying Proposition 15.51 and substituting the values of all parameters c^ we 
obtain the desired formula for correlation functions R n . □ 



6 Bulk limits 

In this section we aim to compute so-called "bulk limits" of the correlation 
functions introduced in Section [5] 

Note that while in the previous sections parameters N, T were always the 
same, through this section they will change. 

We are interested in the following limit regime: Let us fix positive numbers 
So, T, N, t, x. Introduce a small parameter s <C 1, and set 

So = S Q e- 1 +o(e- 1 ), T = fe" 1 + o(e _1 ), N = Ne' 1 + o(e _1 ). 

Consider also integer valued functions ti = t»(e) and x% = Xi(e), i = 1, . . . , n, 
such that 

lim sti(e) = t, lim exi{e) = x, i = 1, . . . , n, 

and pairwise differences U — tj, and Xi — Xj do not depend on s. 
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Then correlation functions R n defined in Theorem 15. fl tend to a limit R n 
which depends on the parameters of the limit regime 5*0, T, N, t, x. 

We consider the region where the limit correlation functions are nontrivial. 
This region is commonly referred to as "bulk" and it is characterized (see e.g. 
|Gor] ) . by the fact that the expression 

-N{N + f ) + {-x + S + N){i+N -x) + £{f + x - S - i) 
2^Jx{-x + S + N)(t + N - x){x + f-S Q -i) 

is strictly between —1 and 1. Denote by <j> = cj)(So,T,N,t,x) the arccosine of 
this expression. 

Let us denote by Qi the triplet (ri,ti,Xi). 

Theorem 6.1. Let n < r2 < ■ ■ ■ < r„, t\ > t-x > • • • > t n and R n (Qi, . . . , Q n ) 
be defined as in Theorem ] 5. 11 Then 

lim R n (Q u ...,Q n ) = det[^ lk ] i)i=1 ,..., n , 
where for (i,j) such that Ti < rj, or Ti = r,j, ti > tj 

1 / a +<**)«•-*. ft dz 

'1- k=n+l 

and for (i,j) such that Ti > rj, ti < tj 



K T=hj _ {i+ciz) "~"- n a +«*-') 



hi <—)"-"-, ft (i 



Tv^bulk 



dz 



/ x{S + N - x) 



(T -t-S +x)(t + N -x) 



J x(i+N-x) 



(T-t-So + x)(S a + N~x) 



Here j± are contours in C joining e l <P( s o,T,N,t,x) an( ^ e i<j>{S ,T,Nd,x) an( ^ 
crossing R± , respectively. 

Comments. The limiting correlation functions 

Rn(Qu ■ • ■ , Qn) = det[K^% J=1 ,...,„ 

are correlation functions of the limit process defined on a fixed admissible sec- 
tion. The proof of the existence of the limit process can be found for instance in 
|Bor[ Lemma 4.1]. The limit process satisfies certain Gibbs property, see |BS| . 
The case r% = r% = • • • = r n = was thoroughly studied in |Gor] . The limit 
process for this case first appeared in [OR] . Our argument is based on the facts 
proved in [Gorj . 
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Proof. Our first goal is to find the limit of the correlation kernels 
lixa K(n,ti(e),Xi(s);rj,tj(e),Xj(e)). 

s— »0 

Let us introduce six auxiliary families of Z x Z matrices or, equivalently, 
operators acting in ^(Z). Set 

{N—l 
E *f(*)*f(y), x,y€X s >\ 
k=0 
for other x, y, 

(recall that 'JJV (x) and X s,t depend on TV" and T, although these indices are 
omitted) 

' (a;,y) = <( fc>w 

for other x, y. 

Observe that P — P' = Id^s.t. Define also 



V^' T ' 5 '*(a: ) i/) = ^>o 



for other a;, y, 



E CW*r(aOn' (v) x e x s ^ y e X s '*- 1 , 



for other it, y, 



^!">>y)=^- n 



E 4+(*)** (*)**' (tf) * e X s -*, y e i^+i.*, 



r N,T,S,t/ 

■''■in = \ 

for other x, y, 

. E41(*)^''(*)*r M (y) *€**«,„€**-* 

VgJ (x,y) = < fe>o 



for other a;, y. 

S.t s,t 

v t± > v s± 



These are (trivial) extensions of finite matrices v t ± , v s ± of Section [2.41 
Now let 



S = Sae^ + oie- 1 ), T = Te' 1 + o^- 1 ), N = Ne' 1 + o^- 1 ), 
t = ie^ 1 + o(e _1 ), x = [xe^ 1 ] + x, y = [is -1 ] + f, 

and send £ to 0. All 6 families tend to some limits. Let us denote these limit 
operators by P(x,u), P'(x,v), V t+ (x,v), V t _(x,v), V s+ (x,v) and %_{>t,v), 
respectively. These limit operators depend on the parameters of limit regime 
N,f,S , i, x. 

As we are dealing with linear operators in ^(Z), it is convenient to employ 
the Fourier transform 

F : J a (Z) ^ L 2 (S 1 ) 1 
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where S' 1 = {zgC:|z| = 1}. We denote the images of our operators under 
Fourier transform by FP, FP' and so on. 

Proposition 6.2. When e — > operators P N - T - S ' t and P' ' ' ' strongly con- 
verge to limits P and P' , respectively. 

FP is the operator of multiplication by the characteristic function of the 
right arc of the unit circle contained between the angles —(f) and <fi. 

FP' is the operator of multiplication by the characteristic function of the left 
arc of the unit circle contained between the angles —(f) and <p. 



Proof. See Section 3.2 in [Got] . □ 

AT T *? / 

Proposition 6.3. When e — > operators V t <' ' ' strongly converge to a limit 
Vt+ ■ FVt+ is the operator of multiplication by the function 



'(T-i-So + xjii+N-x) _ 1 J (S Q + N-£)£ 



V (t + N){T + N -t) V (t + N){T + N -t) 

Proof. See Section 3.3 in [Gorj □ 

A7 T 1 *? t 

Proposition 6.4. When e — > operators V t J ' ' strongly converge to a limit 
Vt~ . FVt— is the operator of multiplication by the function 



'(T-i-So + X^i+N-x) J (S + N-S:)S: 



V (t + N)(T + N-t) V (t + N)(T + N-t) 

Proof. Recall that c t L (k) = c t _J_ _ (fc) and c t L (k) is a real number. Conse- 
quently, V t _' ' ' — (V t+ ' ' '~ )*. Therefore, operators V t J ' ' tend to a limit 
Vt- and FVt- = (FVt+)* . Finally, observe that if A is the operator of multipli- 
cation by z _1 acting in ^(■S' 1 ), then A* is the operation of multiplication by z. 
Thus, FVt- is given by the desired formula. □ 

Proposition 6.5. When e — > operators V s 2 ' ' strongly converge to a limit 
Vs+ ■ FVs+ is the operator of multiplication by the function 



! (T-t-S Q + £)(So + N - x) _ x / (t + N-x)x 



V (So + N)(T + N - So) V (So + N)(T + N - So) 

Proof. Results of Section [2] imply that Vgl ' ' ^ V t+ ' ' ' under the S <-> t. 
Perform this involution, then send e to zero and then switch S <-> t again. 
Proposition 16.31 implies the result. □ 

A7 T 1 S t 

Proposition 6.6. When e — > operators VoJ ' ' strongly converge to a limit 
Vs- ■ FVs- is the operator of multiplication by the function 



'(T-i-So+x)(So+N-x) J (t + N~x)5; 



(So + N)(T + N~S ) V (So + N)(T + N - S ) 
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Proof. Same argument as in Proposition 16.51 □ 



Now we proceed to the proof of Theorem 16.11 The correlation kernel 
K(ri, ti,Xi', Tj, tj, Xj) defines aZxZ matrix or, equivalently, an operator acting 
in 1%(I) by 

tri N _ J ^"'' ' X ' '"•?' ■?' ^' X N,T ' V ^ N,T ' 

1 for other x, y. 

First, suppose that n < Tj or r, = rj, U > tj. In this case K,(x,y) can be 
decomposed in the following way: 

K(X y) = p> N ' T ' S ( r i)> t i . yN,T,S(r 3 ),t ] + l _ v N,T,S(r 3 ),t ]+ 2 yN,T,S(rj),t t 

v N,T,S(r t -l),U . v N,T,S(r t -2),U yN.T, S (r t ),U (1Q) 

i>e rj . i>£ rj . _i ' be ri + i \ I 

This relation readily follows from the definition of K,(x,y) and orthogonality 
relations on functions \& ; ' (x). 

Observe that norms of all factors in (JTUJ) are bounded by 1. ( P' ' ' 2 

are ortoprojection operators; norms of V's are bounded since constants c used 
in their definition are bounded by 1). Consequently, convergence of each factor 
in (fTOf as e — > implies strong convergence of K.(x,y) to the limit operator 

£ = £'($_)*«-*' • V S e Tj ■ V S e rj ^ ■ ■ ■ V S e ri+1 . 

(Indeed, multiplication of operators is jointly continuous on bounded sets in 
strong operator topology.) The image of /C under the Fourier transform is given 
by 

FK = FP'tft)''-' 3 • J] FV Stk+1 . 
Performing the inverse Fourier transform we obtain the formula 
]xmK(ri,ti(e),Xi(e)',rj,tj(e),Xj(e)) = K(xi,Xj) 



1 T l l{T-i-S + x)(i+N-x) J (S + N-5;)x 



2iriJ e -n,\y (t + N)(T + N-i) \ (t + N)(T + N -t) 

r i 

n 



x II I ( T -t-So + x)(So + N-x) z _ ek J (t + N-x)x 



, . , v s (So + N)(T + N-S ) V (S + N)(T + N - S ) , 

dz 



x — 



r-r — const/ 3 ■ COUStry 3 % ■ K;„ 



with 



(T-i-S + x)(i+N-x) 
consU = \ 



(t + N){T + N-t) 
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canst, = KT-i-So+^^ + N-xl 

V (So + N)(T + N-S ) 

(The integration is performed over the left side of the unit circle.) 

The case r, > Tj, ti < tj is quite similar. The only difference is that instead 
of the operators V t J ' ' and V s ± ' ' we have to use in some sense inverse 
operators. Actually, these operators are not invertible and limit operators Vs± 
and Vt- might be non-invertible too, but the difficulties can be avoided if we 
restrict all operators on the images of p N ' T ' S ' t and P. Details of this trick can 
be found in Section 3.3 of [Gorj . 
The answer for n > Tj, ti < tj is 

lim K(n,ti(e),Xi(e);rj,tj(e),Xj(s)) 

e—>0 




1 T / l(T-i-S + x){t + N -x) J [S +N-x)x 



2iriJ e -^\V (t + N)(T + N - 1) V (i + N)(T + N - i) 



-A- / (T-i-S + x){S + N-x) | z _, k J (i+N-x)x 



X 

k=r,+l 



(S + N)(T + N-So) V (So + N)(T + N - S ) 



dz 



= const'r' 3 ■ consQ~ Tl ■ K hulk 



z Xi-Xj + l ^"'-""1 ^.^u 2 iiy , 

where consti and cunsti are as above. (The integration is performed over the 
right side of the unit circle.) 

Since in the determinant for correlation functions the prefactors 
consti 3 const 2 ' cancel out, the proof is complete. 

□ 
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